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We consider two versions of quantum Regge calculus. The Standard Regge 
Calculus where the quadratic link lengths of the simplicial manifold vary 
■ continuously and the Z2-Regge Model where they are restricted to two possible 
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values. The goal is to determine whether the computationally more easily 



accessible Z2 model still retains the universal characteristics of standard Regge 

theory in two dimensions. In order to compare observables such as average 

curvature or Liouville field susceptibility, we use in both models the same 

functional integration measure, which is chosen to render the ^-Regge Model 

particularly simple. Expectation values are computed numerically and agree 

qualitatively for positive bare couplings. The phase transition within the 

^-Regge Model is analyzed by mean-field theory. 
PACS: 04.60.Nc 
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I. INTRODUCTION 



Standard Regge Calculus (SRC) |T| provides an interesting method to explore quantum 
gravity in a non-perturbative fashion ||. The infinite degrees of freedom of Riemannian 
manifolds are reduced by discretization, that is, SRC deals with piecewise linear spaces de- 
scribed by a finite number of parameters. A manifold is approximated by a simplicial lattice 
with fixed coordination numbers, as opposed to the dynamical triangulated random surface 
(DTRS) method || where the coordination numbers are treated as the dynamical degrees of 
freedom. This leaves in SRC the quadratic link lengths q as gravitational degrees of freedom 
which are constraint by triangle inequalities. Since analytical treatments have proven to be 
difficult , this approach has been extensively studied through numerical computer simulations 
during the last ten years ||||. Although the computer codes can be efficiently vectorized, 
large scale simulations are still a very time demanding enterprise. One therefore seeks for 
suitable approximations which will simplify the SRC and yet retain most of its universal 
features. 

The Z 2 -Regge Model (Z 2 RM) || could be such a desired simplification. Here the 
quadratic link lengths q of the simplicial complexes are restricted to take on only the two 
values qi = 1 + ea h < e < e max ,cri = ±1, in close analogy to the ancestor of all lattice 
models, the Lenz-Ising model. To test whether this simpler model is in a reasonable sense 
still similar to SRC, i.e. shares the same universal properties, we study both models in two 
dimensions and compare a number of observables for one particular lattice size. Moreover 
we estimate in both models the critical exponent r\§ of the Liouville field susceptibility by 
performing a finite-size scaling analysis on moderately sized lattices. 

Although some models for 2d-quantum gravity have been exactly solved via the matrix 
model approach J7| and with the help of conformal field theory ||, the relation of those 
approaches to SRC is not yet understood. There are even severe discrepancies between the 
alternative discrete approach, the DTRS method ||, and SRC. Especially the functional 
integration measure in SRC is under heavy debate ||. In an effort to clarify the role of the 
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measure the conventional definition of diffeomorphisms has been employed, assuming that a 
piecewise linear space, i.e. a Regge surface, is exactly invariant under the action of the full 



diffeomorphism group (K| . After a conformal gauge fixing was performed in the continuum 
formalism, it was shown that the evaluation of the non-local Faddeev-Popov determinant 
by using such a Regge regularization leads to the usual Liouville field theory results in the 
continuum limit. All that is based on a description of piecewise linear manifolds with deficit 
angles, not edge lengths, and is mostly taken as an argument that the correct measure of 
Standard Regge Calculus has to be non-local. However, to our knowledge it is not obvious 
that this argument carries over to a discretized Lagrangian, which is formulated in terms 
of fluctuating edge lengths, obeying triangle equalities, and which is not invariant under 
the diffeomorphism group due to the presence of curvature defects: different assignments of 
edge lengths correspond to different physical geometries, and as a consequence there are no 
gauge degrees of freedom in Standard Regge Calculus, apart from special geometries (like 
flat space) [11]. Therefore we do not include a gauge fixing term and rely in this work still 
on a local measure. Our main goal in this present investigation is not to resolve the measure 
question but to explore the phase behavior of the Z 2 RM and its relation to SRC. We will 
show that the discretized Z 2 RM does not suffer from unphysical gauge degrees of freedom. 
If both SRC and Z2RM for certain local measures lie in the same universality class one can 
hope to learn about physical observables using this simplified approach. 

The rest of the paper is organized as follows: In Chapter |T| we briefly review the Standard 
Regge Model as well as the Z 2 -Regge Model. In Chapter |T| we introduce the observables 
and discuss important scaling relations. The details of the Monte Carlo simulations and the 
results are presented in Chapter [IV]. Chapter [V] deals with a mean-field approach for the 
Z 2 RM to discuss the observed phase transition in that specific model. Finally Chapter |VT| 
ends with our conclusion. 
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II. MODELS 



Starting point for both Standard Regge Calculus and the Z 2 -Regge Model is Regge's 
discrete description of General Relativity in which space-time is represented by a piecewise 
flat, simplicial manifold: the Regge skeleton The beauty of this procedure is that it 

works for any space-time dimension d and for metrics of arbitrary signature. The Einstein- 
Hilbert action translates into 

I{q) = A £ V(s d ) - 2(3 ^ 5(s d - 2 )V(s d - 2 ) , (2.1) 

s d s d-2 

with the quadratic edge lengths q describing the dynamics of the lattice, A being the cosmo- 
logical constant, and (3 the bare Planck mass squared. The first sum runs over all (i-simplices 
s d of the simplicial complex and V(s d ) is the ci-volume of the indicated simplex. The second 
term represents the curvature of the lattice, that is concentrated on the (d — 2)-simplices 
leading to deficit angles 5(s d ~ 2 ), and is proportional to the integral over the curvature scalar 
in the classical Einstein-Hilbert action of the continuum theory. The connectivity of the 
edges, in simplicial terminology called the incidence matrix, is fixed from the beginning 
through the simplicial decomposition of the manifold under consideration. Any smooth 
manifold can be approximated by a Regge skeleton with arbitrarily small deficits simply by 
using a sufficient number of links and arranging them appropriately. 

A. Standard Regge Calculus 

In two dimensions Regge's discretization procedure is easily illustrated by choosing a tri- 
angulation of the considered surface. Each triangle then represents a part of a piecewise lin- 
ear manifold. The net of triangles itself is a two-geometry, with singular (non-differentiable) 
points located at the vertices of the net. In the presence of curvature a vector that is parallel 
transported around a vertex experiences a rotation by the deficit angle 5i = 2n — J2tDi ^*(0> 
where 9i(t) are the dihedral angles of the triangles t attached to vertex i. The integral of the 
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scalar curvature over the simplicial complex K in two dimensions is a topological invariant 
due to the Gauss-Bonnet theorem. Its simplicial analogue reads as 

J2 5, = 27r X (M) , (2.2) 

iDK 

where x — 2(1 — g) is the Euler characteristic of the manifold A4 expressed by the number 
of handles g in A4. 

In the exceptional case of flat skeletons one can move a vertex on the surface, keeping all 
the neighbors fixed, without violating the triangle inequalities, such that different configu- 
rations triangulate the same (flat) geometry. This transformation has two parameters and 
is an exact invariance of the action, but does not exist in general. When space is curved 
the invariance is only an approximate one. In the limit of increasing number of links local 
gauge invariance, that is the continuum diffeomorphism group should be recovered. 

A quantization of the above action ( |2.1|) proceeds by evaluating the path integral 



Z = J D[q]e- I( - q) . (2.3) 

In principle the functional integration should extend over all metrics on all possible topolo- 
gies, but, as is usually done, we restrict ourselves to one specific topology, the torus, M. = T 2 . 
Consequently the Euler characteristic x{T 2 ) vanishes in (|2.2|) and the action (|2.1| ) consists 
only of a cosmo logical constant A times the sum over all triangle areas A t . The path-integral 
approach suffers from a non-uniqueness of the integration measure, even the need for a non- 
local measure is advocated. However, some of the proposed non-local measures do not agree 
with their continuum counterparts in the weak field limit, which is a necessary condition 
for an acceptable discrete measure |J. This property however is fulfilled for the standard 
simplicial measure |L2 

/£>[?] = rL /%^(?0, (2-4) 

J HI 

with m G M permitting to investigate a 1-parameter family of measures. The function 
Fniqi) constrains the integration to those Euclidean configurations of link lengths which 



do not violate the triangle inequalities. The positive parameter r/ modifies the triangle 
inequalities to Z 3 < (Zi + Z 2 ) (1 — rj) and Z 3 > \l\ — l 2 \(l + rj), so that very thin triangles are 
suppressed. This is not necessary on theoretical grounds, but will be useful for the Monte 
Carlo evaluation of the path integral. 

Hence the model considered here is characterized by the partition function 

^({*})e- AE ** , (2-5) 

where N± is the number of links and Aj = J2tDiAt/3 denotes the barycentric area with A t 
being the area of a triangle t. A specific choice of the value of the parameter m will be 
discussed in the next section. 

As discussed above, the numerical computations of (|2.5| ) (for technical details see Chapter 
[TVD do not run into the diffeomorphism problem by summing over distinct simplicial lattices 
without fixing a gauge. Still the question arises whether one double-counts some classes of 
geometries in this way, e.g., one may argue from the above that flat geometries are over- 
represented, though simulations give no hints on that. 
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B. Z 2 -Regge Model 

In the Z 2 -Regge Model || the squared edge lengths q\ are allowed to take on only the 
two values 

qi = 1 + tai , < e < e max , o\ = ±1 , (2.6) 

where the parameter e is chosen such that the Euclidean triangle inequalities are fulfilled 
for all qi's, i.e. T.^ = 1 for all configurations {qi}- There exist 2^ different configurations 
and for finite e and link lengths none of them can be transformed smoothly into each other. 
A further nice attribute of the Z 2 RM is its accordance with lattice perturbation theory. As 
described in |Tj|, ( |2.(j| ) can be viewed as weak-field expansion around flat space, implicitly 



having performed a conformal gauge fixing. On each triangle the metric tensor assumes 
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the form ^(A) = (1 + t)5^ v . However a finite e inhibits local conformal transformations. 
Since triangles share links, rescaling lengths on one particular triangle would necessitate 
the same rescaling for the neighboring triangle and so on. What remains only is a global 
transformation. In the continuum limit the local conformal degree of freedom should be 
restored. 

Using Q2.6|) the measure ( ^4|) can be replaced by 

exp[-m^ln(l + eai)\ = exp[-iVim (e) - m i(^) cr i] , ( 2 -7) 

<7;=±1 I CT;=±1 I 

where m = — |me 2 + 0(e 4 ) and mi = m[e + ^e 3 + 0(e 5 )] = mM with 



oo 2i-l i i , 

M = Y. J — ='ln 



f t 2i-l 2 1 - e 

The area of a single triangle t with squared edge lengths qi, q 2 , qi can be expressed as 



2.8) 



.4, 



1 



\{q\ + q2-qi) q2 



13 1 1 3 i 

2^4 + 2^ + a2 + ai ^ e + 2^ ai(r2 + ° xai + <72<Jl ~ ' 



(2.9) 



Expanding A t the series consists only of terms up to a 3 since of = 1. Therefore A t can be 
written as 

A = c (e) + Ci(e)Oi + cr 2 + a;) + c 2 {e){a 1 a 2 + o\Oi + a 2 a{) + c z {e)a x a 2 ai . (2.10) 



Computing the four possible values for the triangle areas and comparing with (|2.10 ) results 
in exact solutions for the coefficients c,-, 
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c = 


32 
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Cl = 




32 




1 


c 2 = 


32 




1 


c 3 = 


32 



2V3 + 3J(l-e)(3 + 5e) + 3 J(l + e)(3 - 5e 



2V3e + J(l - e)(3 + he) - J (I + e)(3 - he) 



2V3 - J(l - e)(3 + 5e) - J(l + e)(3 - 5e 



2^ - 3J(1 - e)(3 + 5e) + 3^(1 + e)(3 - 5e) 



(2.11) 
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Obviously one must have e < | = e max for the triangle areas to be real and positive. In the 
simulations described below we used e = 0.5 but the special choice of e does not affect our 
numerical results except for the value of the critical coupling (A c ) as one can infer from the 
next equation. It is obtained by inserting ( |2.7| ) and ( |2.10| ) into the partition function (|2.5| ) 



Z = Jexp{-J2(2Xc 1 + m^ai - \Y^[c2(cricr2 + c^icri + cr 2 cri) + c^a^ai]} , (2.12) 

ai=±l I t 

with J = exp(— AiV^co — Niitlq) and N 2 denoting the total number of triangles. If we 
view Ui as a spin variable and identify the corresponding link / of the triangulation with a 
lattice site, then Z reads as the partition function of a spin system with two- and three-spin 
nearest-neighbor interactions on a Kagome lattice. A particularly simple form of Eq. (|2.12|) 
is obtained if one chooses m x = — 2Aci, and therefore 

which is henceforth used for the measure in the Z 2 RM as well as in SRC. 



III. OBSERVABLES AND SCALING RELATIONS 

To compare both models we examined the quadratic link lengths and the area fluctuations 
on the simplicial lattice. Additionally we consider the squared curvature defined by 

i? 2 = Ef- (3-1) 

Furthermore the Liouville mode is of special interest because it represents the only degree 
of freedom of pure 2d gravity. The discrete analogue 0, of the continuum Liouville field 



<p(x) = In J g(x) is defined by 

<& = hiA, A = \j2 A t' (3-2) 

6 tDi 

where Ai is the area element of site i. The Liouville susceptibility is then defined by 

X, = (A)m - <0> 2 ] , (3.3) 



where = an d A is the total area |j) . ^From Ref . || it is known that for fixed total 

area A and the scale invariant measure dq/q the susceptibility scales according to 

ln X0 (L) i ~°°c+(2-^)lnL, (3.4) 

with L = \/~A and the Liouville field critical exponent 77^ = 0. 

One can easily derive a scaling relation for the SRC from the partition function (p.5|) 
HI . Rescaling all quadratic link lengths of Eq. (|2]5|) by a factor (, i.e. q' = C -1 ?, yields 



, Jo 



FviWiVe-^W • (3.5) 



Since this is only a change of the integration variables, Z cannot depend on (. Hence we 
have 

dlnZ 



1 = = N 1 (l-m)-\(£A i ) (3.6) 



and find that the expectation value of the area A is fixed to 

which is a useful identity to check the correctness of the simulation results. 

For positive A the scaling relation makes always sense, and we expect the partition 
function ( |2.5| ) to behave well. However, there is also a range for negative values of 
A < \ c ,scai = — M/2ci, where the area expectation value remains positive. For our choice 
of e = 0.5 one obtains M = 0.5493 . . . and c\ = 0.0788 . . ., such that this value equates to 
Kscai = —3.4817. For values of \ csca i < A < the area expectation value is negative and 



hence the partition function (2.5) is ill-defined 



IV. SIMULATION DETAILS AND RESULTS 



We studied the partition functions ( ^.5[ ) and Q2.12p on toroidal lattices with iVo = 16 x 



16 = 256 sites, resulting in N\ = 768 and N 2 = 512. The measure was chosen as a 
function of A according to m — — ^jg 1 = \/\ c ,scai with e = 0.5, i.e., m w — 0.3A. In the 



SRC formulation we used a standard Metropolis algorithm to update every link during a 
lattice sweep. For each A in the range of 1 — 30, starting from an initial configuration with 
equilateral triangles, we first thermalized the system and then collected 200k measurements 
taken every 50 th MC sweep. Error bars were computed by the standard Jackknife method 
on the basis of 10 blocks. The integrated autocorrelation time ta for the area was in the 
range of unity, whereas the integrated autocorrelation time r#2 of R 2 was about 100 — 500, 
and t ^ of was about 100 — 1000, depending on the parameter A. 

For the update of the spins in the Z2RM formulation we also used a Metropolis algorithm. 
Here we performed 100k MC update sweeps measuring every 10 th sweep after an initial 
equilibration period. The integrated autocorrelation times were typically of the order one. 

In our first test runs of SRC with positive A we noticed that our measured area expecta- 
tion value was not in agreement with the scaling relation ( |3.7|) , and that the mismatch was 
growing when we increased A. A way to cure this problem was to implement a new global 
MC move, which we termed the breathing move. It simply consists of rescaling all quadratic 
link lengths q by a factor (, which is the same as scaling the total area by the factor (, 
where ( can be smaller or larger than one. Again, the usual Metropolis criterion was used 
to accept or reject this proposal, and the extent of ( was adjusted to yield an acceptance 
rate of roughly one half. We allowed for ten global breathing moves followed by a complete 
sweep through the lattice with normal link updates. Applying these two moves the area 
expectation value decreases with increasing A in perfect agreement with the scaling relation 
( p.7| ) , see the upper left plot in Fig. |]. For sake of better comparison observables are plotted 
in Fig. [i] and |] relative to the critical cosmological constant, A* = A — A c (= A for SRC). 

With the above global scaling one can also show that the area expectation value, which 
follows from the scaling relation for A < X c , S cah is an unstable equilibrium value. This can 
be seen by considering the ratio of the Boltzmann factors B before and after a proposed 
global move, where we set g new = (q old and A n ™ = (A old , with ( = 1 + k and m — 1 — a. 
Using these abbreviations we get 
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£? ncw _ exp(-ACA old + a £ In C<f ld ) 
^oid _ exp(-AA old + a£lng old ) 

= exp(-Afi:A old + a^ln(l + fi:)) (4.1) 

« exp(- K {AA old -aNx}- + ...). 

The moves are always accepted in the Monte Carlo simulations if > 1. The first term in 
the argument of the exponential is linear in k and constitutes essentially the scaling relation 
( P.7D - The second term quadratic in k is always positive for negative A < A CiSca ;, i.e., for 
negative a. Therefore this term will initially drive the system away from its equilibrium 
value. In the next steps, the linear term amplifies the instability and, depending on the 
arbitrary initial sign of k, the area will tend to zero or infinity. 

However, even with the refined update scheme we observed that the SRC system ther- 
malizes extremely slowly for very small A and therefore display only statistically reliable data 
points for A > 1 in the plots on the l.h.s. of Figs. [l| and 0. One also expects that expectation 
values of the average squared link length q will decrease and that expectation values of the 
squared curvature will increase with growing A which is indeed confirmed by the simulation 
results displayed in Fig. [TJ. The Liouville field (0) and the associated susceptibility Q3.3| ) are 
shown on the l.h.s. of Fig. The peak in the Liouville susceptibility can be expected to 
move towards A = for L — > oo due to the diverging area at this point. 

The corresponding quantities of the Z 2 RM are shown for comparison on the r.h.s. of 
Figs. [I] and 0. Since it is a crude approximation a quantitative agreement with SRC cannot 
be expected. We see, however, that the two formulations yield for most quantities the 
same qualitative behavior for positive couplings A*. Whereas the SRC becomes ill-defined 
for negative A*, the Z 2 RM as an effective spin system is well-defined for all values of the 
cosmological constant. In all quantities of the Z 2 RM we observe the signature of a crossover 
or phase transition, which is particularly pronounced in the Liouville susceptibility. From 
the peak location we read off a critical coupling A c ~ —11. The phase transition in the 
Z 2 RM might be viewed as some relic of the transition from a well- to an ill-defined regime 
of SRC. To get some idea on the nature of the phase transition we looked at the histograms 
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of the total area and the squared curvature in Fig. |3|. They show a single-peak structure at 
A c which would hint at a continuous transition. 

By performing several runs, we further have monitored the scaling of the Liouville sus- 
ceptibility maxima in the Z2RM for lattice sizes up to L = 256. The double-logarithmic plot 
in Fig. [| shows that the scaling behavior fl3.4| ) is governed by a critical exponent 77^ ~ 2. 
We have performed a finite-size scaling analysis of the Liouville susceptibility on lattice sizes 
L = 4, 8, 12, 16, 24, and 32 also for SRC. Because of the transition to an ill-defined phase 
making computations in its vicinity difficult we employed in all simulations a fixed value 
A = 1 closest to the critical coupling (see Figs. [I] and 0). The Liouville-field exponent gives 
the behavior of the Liouville field in a specific phase, here the well-defined phase, and hence 
should be independent of A. After an equilibration period we performed 20k - 100k measure- 
ments of the Liouville field and computed the associated susceptibility ( |3.3|) . Our results 
presented in Fig. f| indeed show that we obtain the same value of 77^ ~ 2, consistent with 
the results found in the Z 2 RM. Furthermore, first results demonstrate scaling of physically 
relevant quantities with the number of allowed values for o\ indicating universality between 



the ZoRM, SRC and models "in between" HI 



Thus we observe strong disagreement with the Liouville value 77^ = 0, which has been 
found for SRC with the scale invariant measure dq/q and fixed area constraint f|. There 
the same lattice transcription as in the present study was employed and it seems reasonable 
to attribute the value 77^ = to the scale invariance of the measure, whereas our measure 
dq/q m with m given in ( |2.13j ) explicitly breaks scale invariance. One therefore might not be 



surprised to observe a different value of r\^. 

V. MEAN-FIELD CALCULATION FOR THE Z 2 -REGGE MODEL 

The variational derivation of mean-field theory starts with a partition function 

z = E»p[-WW)] • (s- 1 ) 
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The comparison with equation (|2.12|) leads to 



Z = ^exp -f3{- [Qi + (oi + o 2 + 0-3 + 0-4) + Q 3 (<7i<t 2 + 0-30-4) <ti} , (5.2 
M V l ) 



with/3 = -|Ac 2 ,ai = ^+ AC2 



and Q3 = The spins cr,-, j = 1, . . . , 4, correspond to 



the squared edge lengths qj = 1 + eaj of the two triangles which have qi and hence 0% in 
common: 




Now the Hamiltonian is divided into two parts, H = Hq + H\, where the choice of Hq is 
only governed by the requirement that it should be possible to evaluate the corresponding 
partition function Z analytically. With the mean-field Ansatz 



H = - T X*i 



(5.3) 



we get 



Z = £ exp f M ]T °) = (e Mo + e 

{a} V i J 



'•) V: 2" l cosh" 1 (M ) (5.4) 
for the partition function. According to Peierl's inequality we obtain for 



Z = e~ pF = Z (e-^ H - H ^) > Z e-^ H - Ho ^ = 



(5.5) 



a lower bound with the mean-field free energy 



- pFmf = In Z - P(H - H ) . (5.6) 
The term (. . .) represents the average with respect to the system described by H , i.e. 
(H - H ) = -— 2 S [Qi + (°"i + o- 2 + 0-3 + cr 4 ) + Q 3 (vi(T2 + 0-30-4) - 



Z 

Mo" 

' P - 



(5.7) 
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Because the thermal average or "mean" value of the spin results in 



Ho = e M + e -M = tanh ( M 0) , (5-8) 

the different terms in ( |5.7| ) give 

(5>,) = N 1 tanh(M ) , 
l 

<£[(<n + ^2 + cr 3 + ^H) = 4iV 1 (tanh(M )) 2 , (5.9) 
l 

(£[(^2 + ^4)^]> = 2iV 1 (tanh(M )) 3 . 
i 

Thus we get a /3- and Mo-dependent expression for the mean-field free energy per link, 
(3f mf = ^ = - In 2 - In (cosh(M )) - 4/3(tanh(M )) 2 - 2Q 3 /5(tanh(M )) 3 + 

iVx 

+(M - /30i)tanh(M Q ) , (5.10) 



that has to be minimized (Z > exp[— /3F m f]). Differentiating ( |5.10 ) with respect to Mq leads 
to the mean-field equation 

tanh(M ) + ^(tanh(M )) 2 = ^ - f (5.11) 

for the free energy. We restrict ourselves again to a calculation without external field Q\. 



A trivial solution of (|5.11| ) is Mq = yielding (3f m f = — In 2, which is the stable solution 
for small (3 (high "temperature"). With increasing (3 a second minimum develops whose 
free energy is eventually lower than that of the M Q = solution. This corresponds to the 



first-order phase transition expected due to the cubic term in ( |5.10|) . To locate the transition 



point we proceed as follows. For a given M value (< 0) we can read off f3 directly from 
Eq. ( |5.11|) , without solving any non-linear equation. By inserting Mq and (3 in ( |5.10|) we 
obtain immediately the free energy which can then be plotted versus (3 as in Fig. |^. This 
curve is double valued because a given M value can correspond to the minimum (lower 
branch) or the maximum (upper branch) of (3f m f at a fixed 13. It should be remarked that 
for (3 > 0.125 there is another meta-stable solution with M > which is not displayed here. 
With decreasing (3 the free energy of the non-trivial solution increases until at (3 c ,mf it hits the 
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value of the M = solution. For (3 < fl c ,mf the non-trivial solution first becomes meta-stable 
and then disappears completely (at the point where the minimum- and maximum-branch 
merge). In this way it is straightforward to extract the transition point as 

Pc,mf ~ 0.1174 corresponding to X c ,mf ~ —8.0 with e = 0.5 . (5-12) 

At this value the free energy shows a structure with two minima of identical height (cf. 
Fig. P), which reconfirms that the mean-field calculation for the Z 2 -Regge Model predicts a 
(weak) first-order phase transition. The value ( |5.12| ) is slightly below the numerical estimate 
of A c ~ —11; however, it is a well-known property of mean- field theory to underestimate 
the critical coupling. In fact, the ratio between the mean-field and the numerical value, 
Xc : mf « 0.73, is of the same order of magnitude as the corresponding quantity for the 
exactly solvable two dimensional Ising model on a square lattice where ( A Ac,mf t )j ~ 0.567. 

VI. CONCLUSIONS 

The aim of this work was to examine whether the Z 2 RM allowing for two discrete edge 
lengths is an appropriate simplification of the SRC with continuously varying edge lengths 
in two dimensions. We studied both models by means of Monte Carlo simulations and found 
that the simpler Z 2 RM qualitatively reproduces the behavior of physical observables like the 
Liouville field or the squared curvature for the bare coupling A > 0. A finite-size scaling 
analysis of the Liouville susceptibility yields for both models the same critical exponent 
rj^ m 2. This hints at a continuous phase transition which is also necessary to perform a 
continuum limit. Concerning universality it is a non-trivial result, because in the SRC model 
with scale invariant measure a value of = was found. 

To obtain further insight, mean-field theory was applied to the Z 2 RM indicating a weak 
first-order phase transition. This would prevent one to gain important information about the 
continuum theory. However, the numerical simulations indicate that fluctuations, which are 
neglected in mean-field theory, soften the true nature of the phase transition to second order. 
The details of this transition such as critical exponents etc. are still left to be determined. 
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An interesting question is the influence of allowing for more than two link lengths and the 
convergence of the properties of such an extended Z 2 RM to those of SRC in two dimensions 



TM. With additional degrees of freedom the situation might resemble the more involved 



four- dimensional case where one has to deal with 10 edges per simplex and the non-trivial 
Einstein-Hilbert action J2tDi^t^t with 50 triangles t per vertex % in Eq. ( |2.1| ). Thus the 
action I(q) takes on a large variety of values already for Z2RM and therefore SRC can be 



approximated more accurately 16 



ACKNOWLEDGMENTS 

A. H. acknowledges support from an Erasmus grant during his stay at the FU Berlin, 
where this project was started. J. R. was supported by Fonds zur Forderung der 
wissenschaftlichen Forschung under contract P11141-PHY. W. J. thanks the Deutsche 
Forschungsgemeinschaft for a Heisenberg Fellowship, and also acknowledges partial sup- 
port by the German-Israel- Foundation (GIF) under contract No. 1-0438-145.07/95. Parts 
of the numerical simulations were performed on the North German Vector Cluster (NVV) 
under grant bvpfOl. 



16 



REFERENCES 



[1] T. Regge, Nuovo Cimento 19, 558 (1961). 

[2] N.H. Christ, R. Friedberg, and T.D. Lee, Nucl. Phys. B202, 89 (1982); J. Cheeger, 
W. Miiller, and R. Schrader, Comm. Math. Phys. 92, 405 (1984); H.W. Hamber, Nucl. 
Phys. B (Proc. Suppl.) 99A, 1 (1991); J. Ambj0rn, in Proceedings of the 1994 Les 
Houches Summer School, Session LXII; J. Ambj0rn, Nucl. Phys. B (Proc. Suppl.) 42, 3 
(1995); S. Catterall, Nucl. Phys. B (Proc. Suppl.) 47, 59 (1996); D.A. Johnston, Nucl. 
Phys. B (Proc. Suppl.) 53, 43 (1997). 

[3] V. Kazakov, JETP Lett. 44, 133 (1986); Phys. Lett. A119, 140 (1986). 

[4] H. Hamber, Nucl. Phys. B400, 347 (1993); H. Hamber and R. Williams, Nucl. Phys. 
B435, 361 (1995); W. Beirl, E. Gerstenmayer, H. Markum, and J. Riedler, Phys. Rev. 
D49, 5231 (1994); W. Beirl, H. Markum, and J. Riedler, Phys. Lett. B341, 12 (1994); 
B. Berg and B. Krishnan, Phys. Lett. B318, 59 (1993); W. Beirl, B. Berg, B. Krishnan, 
H. Markum, and J. Riedler, Phys. Lett. B348, 355 (1995); C. Holm and W. Janke, Nucl. 
Phys. B477, 465 (1996); Phys. Lett. B375, 9 (1996); Phys. Lett. B390, 59 (1997). 

[5] M. Gross and H. Hamber, Nucl. Phys. B364, 703 (1991); C. Holm and W. Janke, Phys. 
Lett. B335, 143 (1994). 

[6] T. Fleming, M. Gross, and R. Renken, Phys. Rev. D50, 7363 (1994); W. Beirl, H. 
Markum, and J. Riedler, Int. J. Mod. Phys. C5, 359 (1994); W. Beirl, P. Homolka, B. 
Krishnan, H. Markum, and J. Riedler, Nucl. Phys. B (Proc. Suppl.) 42, 710 (1995). 

[7] V. Kazakov, Phys. Lett. B150, 282 (1985); F. David, Nucl. Phys. B257, 45 (1985); J. 
Ambj0rn, B. Durhuus, and J. Frohlich, Nucl. Phys. B257, 433 (1985); V. Kazakov, I. 
Kostov, and A. Migdal, Phys. Lett. B157, 295 (1985). 

[8] V. Knizhnik, A. Polyakov, and A. Zamolodchikov, Mod. Phys. Lett. A3, 819 (1988). 



17 



J. Ambj0rn, J. Nielsen, J. Rolf, and G. Savvidy, Class. Quant. Grav. 14, 3225 (1997); 
H. Hamber and R. Williams, |hep-tli/9708'0T9| . 



[10] P. Menotti, Nucl. Phys. B (Proc. Suppl.) 63, 760 (1998); P. Menotti and P. Peirano, 
Phys. Lett. B353, 444 (1995); Nucl. Phys. B473, 426 (1996); Nucl. Phys. B (Proc. 
Suppl.) 53, 780 (1997); Nucl. Phys. B (Proc. Suppl.) 57, 82 (1997); Nucl. Phys. B488, 
719 (1997). 

[11] J.B. Hartle, J. Math. Phys. 26, 804 (1985); J. Math. Phys. 27, 287 (1986). 

[12] H. Hamber and R. Williams, Nucl. Phys. B248, 392 (1984); Nucl. Phys. B269, 712 
(1986). 

[13] A. Jevicki and M. Ninomiya, Phys. Lett. 150B, 115 (1985); M. Rocek and R. Williams, 
Phys. Lett. B104, 31 (1981); Z. Phys. C21, 371 (1984). 



[14] E. Bittner, H. Markum, and J. Riedler, Nucl. Phys. B (Proc. Suppl.), in print; [hcp- 
| lat/9809155. 



[15] J. Binney, N. Dowrick, A. Fisher, and M. Newman, The Theory of Critical Phenomena 
(Oxford University Press, Oxford, 1992). 

[16] W. Beirl, A. Hauke, P. Homolka, B. Krishnan, H. Kroger, H. Markum, and J. Riedler, 
Nucl. Phys. B (Proc. Suppl.) 47, 625 (1996); P. Homolka, W. Beirl, P. Homolka, H. 
Markum, J. Riedler, and H. Kroger, Nucl. Phys. B (Proc. Suppl.) 53, 769 (1997). 



18 



FIGURES 



A 
< 

V 



A 
V 



4 
3.5 

3 
2.5 

2 
1.5 

1 

0.5 


1000 
800 
600 
400 
200 


50 p 

40 - 
30 



+ 




- 























10 



X' 



10 



v 20 



10 



+ + 



10 



20 



20 



20 



30 









+ 






+ 






+ + + 








+ 


+ 



30 



30 




FIG. 1. Expectation values of the area A, the average squared link length q, and the squared 
curvature R 2 as a function of the distance to the critical cosmological constant A* for the Standard 
Regge Calculus (left plots) and the Z2-Regge Model (right plots). Nq is the total number of vertices. 
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FIG. 2. Expectation values of the Liouville field (f> and the Liouville field susceptibility x<j> as a 
function of A* for the Standard Regge Calculus (left plots) and the Z2-Regge Model (right plots). 
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FIG. 3. Histograms of the total area and the squared curvature per vertex at the transition 
point of the Z2-Regge Model. 
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FIG. 4. Scaling of the of the Liouville field susceptibility x<p as a function of the lattice size L 
for the SRC (left) at A = 1 and the Z 2 RM (right) for max(x<p)- 





0.005 r 




- 




-0.005 


eg 




c 


-0.01 


+ 




E 

M— 


-0.015 - 








-0.02 




-0.025 




-0.03 - 




0.115 0.12 0.125 

P 



0.13 



FIG. 5. Value of (3f m f + In 2 at the two minima under variation of (3. 
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FIG. 6. Free energy 0f m f + In 2 as a function of Mq at the mean- field transition point. 
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